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On the basis of various Direct Numerical Simulations (DNS) of turbulent channel flows the follow- 
ing picture is proposed, (i) At a height y from either wall, the Taylor microscale A is proportional to 
the average distance £ s between stagnation points of the fluctuating velocity field, i.e. X(y) = B\l a (y) 
with Bi constant, for 8 V <C y < 8, where the wall unit <5„ is defined as the ratio of kinematic viscos- 
ity v to skin friction velocity u T and <5 is the channel's half width. (%%) The number density n a of 
stagnation points varies with height according to n 3 — where y+ = y/& v and C s is constant 

in the range 5 V -C y < S. (iii) In that same range, the kinetic energy dissipation rate per unit 
mass, e, equals where E+ is the total kinetic energy per unit mass normalised by u\ and 

k s = Bi/Cs is the stagnation point von Karman coefficient, (iv) In the limit of exceedingly large 
Reynolds numbers TZe T = 5/S v , large enough for the Reynolds stress —{uv) to equal u\ in the range 
5 V <C y <C S, and assuming that production of turbulent kinetic energy balances dissipation locally 
in that range and limit, the mean velocity U+, normalised by u T , obeys ~ m that same 

range, (v) It follows that the von Karman coefficient k is a meaningful and well-defined coefficient 
and the log-law holds in turbulent channel/pipe flows only if E+ is independent of y+ and !Ze T in 
that range, in which case k ~ n a . (vi) In support of -§^U+ ~ | ^jj, DNS data of turbulent channel 



flows which include the highest currently available values of TZe T are best fitted by E+ ~ §£?4j/. 
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and jj^-C/+ ~ ^ A U+ 1 2//1S with B4 independent of y in 8 V <C y <C S if the significant departure from 
— (uv) ~ u% at these !Ze T values is taken into account. 



I. INTRODUCTION 



The mean flow profile of turbulent boundary layers is 
very widely taken to incorporate an extensive log-law re- 
gion for a very broad range of turbulent wall-bounded 
flows. As a result, the log-law and its so-called von 
Karman constant are a central part of most engineering 
turbulence models. However, renewed interest and new 
measurements over the past ten to fifteen years have led 
to debates (i) on the form of the mean velocity profile 
(is it a log-law or a power-law with very weak power ex- 
ponent?), (ii) on its scalings with Reynolds number and 
(iii) on its dependence on or independence of overall flow 
geometry (see the Theme Issue of Phil. Trans. R. Soc. 
Lond. A (2007), volume 365, on "Scaling and structure 
in high Reynolds number wall bounded flows" ) . Fittings 
of new mean flow data with a log-law lead to a variety 
of values for the von Karman constant k: as low as lie 
for channel flows as high as 0.43 for pipe flows Q 
and k ~ 0.38 for zero-pressure-gradient boundary layer 
flows 0. Due to this non- universality the von Karman 
constant is being renamed von Karman coefficient Q. 

These seemingly small departures from the classically 
accepted Q| k ~ 0.41 value matter because of the reliance 
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that so many turbulence models have on the log-law. In 
a personal communication in early 2009, Phillip Spalart 
mentions that "the recent fluctuations in preferred k val- 
ues leads to changes of the order of 1% in the drag pre- 
dicted by CFD for an airplane with typical RANS mod- 
els" . In terms of fuel economy sustained over years, this 
1% drag change is significant. 

All these are issues which in one way or another would 
have been solved if the turbulence closure problem had 
been solved. In this respect, the main goal of the present 
paper is not to find definitive answers to these specific 
issues, nor is it to find closure relations between different 
order statistics of flow quantities Q . Instead, our main 
goal is to determine relations between mean flow quan- 
tities, such as the mean flow and the dissipation rate of 
kinetic energy, and the underlying flow field topography 
of the fluctuating velocities. Recent results obtained for 
homogeneous isotropic turbulence indicate that it is not 
impossible to find relations between bulk flow statistics 
and the underlying topography of the fluctuating veloc- 
ity field. Mazcllicr and Vassilicos @ related the dissi- 
pation constant C e to the number of zero-crossings of 
velocity fluctuations and were able to take account the 
non-universality of C e by their formula. Goto and Vas- 
silicos went one step further and related C e to the 
number of stagnation points of the velocity fluctuations. 
The stagnation points are the most basic aspect of field 
topography and they have a multi-scale spatial distribu- 
tion [8j. As such, they are a convenient and well-defined 



2 



concept, unlike "eddies". They have also proved useful 
for understanding turbulent pair diffusion [9j. 

In this paper, we present a phenomenology based on 
the underlying topography of the fluctuating velocity 
field which relates the mean flow profile to the multiscale 
structure of stagnation points of the velocity fluctuations. 
We use Direct Numerical Simulations (DNS) of various 
fully developed incompressible turbulent channel flows to 
validate our new approach and propose a resulting new 
starting point for a new intermediate asymptotic analysis 
of the mean flow profile of turbulent channel/pipe flows. 

With the exception of the highest Reynolds number 
DNS channel flow data [To| which we use towards the 
end of this paper, the Reynolds numbers considered here 
range between low to moderate (though, of course, always 
large enough for the flow to be turbulent). In terms of 
lZe T = = 8j8 v where u T is the skin friction velocity, 
8 the channel half-width, v the fluid's kinematic viscosity 
and the wall unit <5„ = v/u r , the highest Reynolds num- 
ber DNS data correspond to TZe T ~ 950 and 2000. 
The Reynolds numbers of our own DNS data range be- 
tween lZe T ~ 110 and 400. This is too low for a direct 
assessment of the log-law but appears to be sufficient for 
the new approach to turbulent mean flow profiles which 
we propose here and which is based on stagnation points 
of the fluctuating velocity field. 

The paper is organised as follows. In section II we de- 
scribe our DNS of turbulent channel flows and in section 
III we present some of the conventional statistics which 
are obtained from our simulations. In section IV we in- 
troduce the stagnation point approach and its applica- 
tion to turbulent channel flows. The phenomenology and 
the mean flow properties implied by the results obtained 
from the application of this approach to our DNS are 
expounded in sections V and VI. Finally, some analysis 
of the highest Reynolds number DNS channel flow data 
fioj is presented in section VII before summarising our 
conclusions in section VIII. 



II. DNS OF TURBULENT CHANNEL FLOW 

We solve the non-dimensionalised incompressible 
Navier-Stokes equations in Cartesian coordinates 

V-tt = 

1 1 (1) 
d t u + - V(u (8 u) + (u ■ V)u] = - Vp + — — Ait 

2 7ce c 

where lZe c = U c 8/v is the Reynolds number based on 
U c = and Ub is the bulk velocity of the flow kept 
constant in time. 

We use the code of Laizet & Lamballais [ll| where spa- 
tial derivatives are estimated using a sixth-order com- 
pact finite-difference scheme and equations (JlJ are nu- 
merically integrated with a frational step method using 
a three-stage third-order Runge-Kutta scheme. The frac- 
tional step method projects the velocity field to a diver- 
gence free velocity field and the Poisson pressure equation 



is solved in Fourier space with a staggered grid for the 
pressure field. The staggered grid for the pressure was 
used for numerical stability purposes as was the skew- 
symmetric implementation of the non-linear term in the 
Navier-Stokes equation (see eqn. (Q])). The grid stretch- 
ing technique maps an equally spaced co-ordinate in the 
computational space to a non-equally spaced co-ordinate 
in the physical space, in order to be able to use Fourier 
transforms in the inhomogeneous wall-normal direction 

To simulate incompressible channel flow turbulence we 
adopted periodic boundary conditions for u = (u,v,w) 
in the x and z directions except at the walls at y = and 
y = 28 where the boundary conditions are either u = 
or borrowed from studies of flow control schemes aimed 
at drag reduction [13, [l4[. The mean flow is in the x 
direction (i.e. (v) = (w) = but (u) ^ 0, where the 
angle brackets (.), in this study, denote averages in the 
x and z homogeneous directions and time expect when, 
in section IIVI stagnation points of u — (u) arc sought, 
in which case the average (it) is only over x and z) and 
the bulk velocity Ub in that direction was kept at the 
same constant value at all times by a control procedure 
which adjusts the mean pressure gradient — d(p)/dx at 
each time step. The choice of Ub is made in accordance 
with Dean's formula lZe T ~ O.lWTZe^ 8 [l5[ for a given 
choice of lZe T . The number of grid points N x , N y and N z 
in the x, y and z directions as well as the domain sizes 
L x and L z in the x and z directions are given in Table 1. 
The domain size in the y direction is of course L y = 28. 

We use different near-wall forcings and boundary con- 
ditions at the walls so as to demonstrate how our stagna- 
tion point approach accounts for the way that different 
wall actuations modify the mean flow profile. Specifically, 
we considered the following three control schemes: 

(i) u = at the walls with forcing f(y) = 
(~Asm(2iry/A)H(A - y),0,0) near the y = wall and 
similar forcing near the y = 28 wall [r| where H is the 
Heaviside function, A = 0.16U^/S ~ u^/S^ and A = 11<5„ 
(case Al). The forcings are applied to the Navier-Stokes 
momentum equations ([1]). This scheme corresponds to 
a steady wall-parallel forcing localised within eleven wall 
units from the walls and uniform in the direction parallel 
to them. This force field averages to zero if integrated 
across the channel; it decelerates the flow closest to the 
wall but accelerates it in the immediately adjacent thin 
region. 

(ii) u — (0,acos(a(x — ct)),0) at the wall [l4| with 
a/U c = 0.05, a/8 = 0.5 and c = -2U C (case A2). This 
boundary condition corresponds to a blowing-suction 
travelling wave on the wall. 

(Hi) u = at the walls and v(x,yd, z,t) replaced by 
— v(x, yd, z, t) at all (x, z) points on the planes yd = 10(5„ 
and yd — 28 — 10(5„ (case A3). This corresponds to a 
computational control scheme whereby the normal ve- 
locity at a distance yd from the walls is made to change 
sign at every time step. 

The numerical parameters of our computations are 



given in Tabic I. Note that the net mass flux through 

TABLE I: (Color online) Parameters for the DNS of 
turbulent channel flow. The term "Forcing" refers to 
wall or near- wall actuations. 



Case 


Forcing 


TZe c 


He-r 


L x 




Nx x N y x N z 


A 


No 


4250 


179 


4ttS 


4x5/3 


200 x 129 x 200 


Al 


Yes a 


4250 


114.4 


4n5 


4tv5/3 


200 x 129 x 200 


A2 


Yes 6 


4250 


222.3 


4tvS 


4tt5/3 


200 x 129 x 200 


A3 


Yes 


4250 


141.6 


4tv5 


471-5/3 


200 x 129 x 200 


B 


No 


2400 


109.5 


4tv5 


2tyS 


100 x 65 x 100 


C 


No 


10400 


392.6 


2tvS 


71 5 


256 x 257 x 256 



"Xu et al. \V& 
b Min et al. QJ 

the wall is zero in all the cases considered here and that 
the y-integrated momentum balance 

(uv) = u T (l - -) 



(2) 



holds for all y in all cases except with Al forcing where 
it holds for A < y < 25 — A. All our non- forced compu- 
tations have been validated against previously published 
databases SUS- 



HI. CONVENTIONAL DNS RESULTS 

When lZe T 3> 1 one might expect an intermediate re- 
gion 5 V -C y <C 5 where production balances dissipation 
locally [HJ], i.e. — (uv)-^(u) ~ e. The idea of such an 
intermediate region is supported by our DNS results (see 
Fig. [T]) which suggest that 



B 2 ■■■ 

tends to 1 as lZe T 



V/e 



(uv) d 



(u)/e 



(3) 



oo in this intermediate region (a 
recent paper [19| proves this asymptotic result by as- 
suming, however, that the mean flow has a logarithmic 
shape in the intermediate region) and that this region 
where this approximate balance holds also increases as 
lZe T increases. The slight discrepancy away from B 2 — 1 
at these moderate Reynolds numbers is well known and 
agrees with other previously published DNS results [IJ . 

In this intermediate region, equation @ implies 
— (uv) — u^. as y/S — >• and y+ = y/S v —5- oo, assuming 



<iii 



U + (where U. 



+ 



(u)/u T ) does not increase 



that 

faster than y+ with p > 1 in this limit. It then follows 
that in this intermediate equilibrium region, 

d(u) u T 



implies 



(4) 



ny Ay 
as TZe T ^> 1 . At finite Reynolds numbers the equation for 
the mean shear in ((4]) should be replaced by ~ J| ^ 
where 



Bo 



(5) 



Case A 
Case A1 
Case A2 
Case A3 
Case B 
Case C 




400 



FIG. 1: (Color online) Profile of the production to 
dissipation ratio. Note the existence of an approximate 
equilibrium layer which grows with lZe T and where 
production approximately balances dissipation. 



Note that even though B 2 and B3 may tend to 1 as 
lZe T 3> 1, they arc definitely different from 1 and even 
functions of y+ and y/5 at finite values of lZe T . 

The mean flow profiles show clear impacts of the con- 
trol schemes on the mean flow (see Fig. For our 
various control schemes at the same lZe c = 4250, the 
skin friction decreases as a result of both case Al and A3 
but increases when the control scheme A2 is applied (see 
Table I). This observation agrees with Fig. [2] where mean 
flow values for cases Al and A3 are higher than for case 
A (no control scheme), and mean flow values are lower 
for case A2 than for case A. 
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FIG. 2: (Color online) Mean velocity profiles. For 
comparison we also plot best log-law fits, o: U+ = (/+. 

14.2,---: U+ = 1 ± s y+ + 0.0, - - -: 
11.2, — : U + = + 5.2 



u + = oky+ 



With reference to the log-law scaling U+ = — log y++B 
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which results from integration of (|4|) if 1/k is independent 
of y, we plot the coefficient y-^U + versus y + in Fig. [3] 

and the coefficient B = E/+— (y-^Uj^j log y + versus y + in 

Fig. g]for all the six different DNS cases of Table I. Note 
that y-§^U + is usually refered to as 1/k but is in fact 
B 2 /(B 3 k) in the present context where k is defined by 
the left-hand side equation in (j4]) . It is only if B2 and B3 
both equal 1 in the equilibrium layer, as may be the case 
when Tle T » 1, that ~ § % yields £ <«) ~ 2* 

and that j/^£/+ becomes 1/k in the equilibrium layer. 

The values of B are affected by the various control 
schemes (see Fig. 2|) in a way consistent with the obser- 
vations made two paragraphs earlier (higher values of B 
for cases Al and A3 than for A and lower for case A2). 
However, it is hard to conclude on the validity of the log- 
law from these results and in particular from the plot in 
Fig. [3] which clearly shows a significant dependence on 
near- wall conditions, lZe T and y+. It may be that the 
log-law is not valid at all or it may be that the log-law 
is not valid unless the Reynolds number is sufficiently 
high, definitely higher than the Reynolds numbers of our 
simulations. 




Case A 
Case A1 
Case A2 
Case A3 
Case B 
Case C 
1/k = 1/0.41 




300 350 400 



FIG. 3: (Color online) The inverse von Karman 
coefficient = Vj^U + versus y + . Taking the definition of 
k to be given by the left-hand equation in (|4]) it is really 
B 2 /(B 3 k) which is plotted against y + . The effects of 
the various near- wall actuations are significant. 



IV. THE STAGNATION POINT APPROACH 

As our direct DNS study of the mean flow equation in 
(@| does not yield clear results, we chose instead to inves- 
tigate the validity of the dissipation equation in (Q| . For 
this we make use of the stagnation point approach which 
has been used recently to show how the number density of 
stagnation points in high Reynolds number Homogeneous 
Isotropic Turbulence (HIT) determines salient properties 




50 100 150 200 250 300 350 400 

FIG. 4: (Color online) B = U+ - (2/3^+) \ogy+ 
versus y + for the six different DNS cases in Table I. 



of turbulent pair diffusion J9| and kinetic energy dissipa- 
tion rate per unit mass ^ n particular, a generalised 
Rice theorem was recently proved Q for high Reynolds 
number HIT which states that the Taylor microscale is 
proportional to the average distance between neighboring 
stagnation points. This average distance is defined as the 
— 1/d power of the number density of stagnation points 
which are points in the d-dimensional space of the flow 
where the turbulent fluctuation velocity is zero. 

The generalised Rice theorem Q for high Reynolds 
number HIT holds under two main assumptions: (i) sta- 
tistical independence between large and small scales and 
(u) absence of small-scale intermittency effects. The 
question which arises in the context of the present work 
is whether it also holds in some region of turbulent chan- 
nel flows. For other cases where measures and concepts 
from HIT impact wall-bounded turbulence please refer to 

MM. 



To obtain some insight into this question by way of 
our DNS, we consider stagnation points of the turbulent 
fluctuation velocity field u' = u — (u), i.e. points where 
all components of the velocity fluctuations around the 
local mean flow are zero. A 3D plot of these points for 
an instant in time in our DNS channel is presented in Fig. 
[5] We use fourth-order Lagrangian interpolation and the 
Newton-Raphson method to locate these points. Details 
on how these points are found are given in the Appendix 

D 

We define N s (y + ) to be the total number of these stag- 
nation points in a thin slab parallel to and at a dis- 
tance y from the channel's y = wall (this number can 
also be defined with respect to the y = 26 wall, but 
here we only consider < y < S without loss of gen- 
erality). The dimensions of this slab are L x x 8 y x L z 



with slab thickness S v 



The average distance be- 



tween stagnation points at a height y from the wall is 
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FIG. 5: (Color online) Points where u' = u — (u) = 
for case C at a given instant in time 
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A Taylor microscale X(y) can be defined 



from e{y) = 2v(sijSij) 



where 



is the fluctuat- 



ing velocity's strain rate tensor and E(y) = ^(\u\ 2 ). The 
question raised is whether a region of turbulent channel 
flow exists for TZe T 3> 1 where 



X(y) = B&iy) 



(6) 



with B\ independent of y and Reynolds number. The 
answer provided by our DNS is that B\ is indeed approx- 
imately constant over an intermediate range 8 V <C y < 8, 
but not perfectly so as our plots in Fig. [5] attest to. It 
is worth noting that this constancy of B\ appears to be 
better defined for cases A, B and C where there is no wall 
or near-wall actuation. Hence, a small discrepancy away 
from B\ = Const may be achieved as a result of those 
different wall-forcings. However, part of the even smaller 
discrepancy in cases A, B and C might be accountable 
to neglected small-scale intermittency effects which, in 
the case of high Reynolds number HIT, are known to 
manifest themselves as a weak Reynolds number depen- 
dence on Bi @. In the case of wall-bounded turbulence, 
small-scale intermittency effects could therefore manifest 
themselves as a weak dependence of B\ on local Reynolds 
number y + = y/S^ (see also [Hj]). However, we leave this 
refinement for future studies. 
Using © and £ s = 



we write 



v2E 
3^ 



v 2E 



2E 



3 B{P S 3 B(L X L Z 



v2E . 
3 B\ 



(7) 



where we have introduced the number density of stag- 
nation points n s = N S /(L X L Z 8 V ). Combining this last 
equation with —(uv)-^{u) = B2£ and using 

as well as C 



dy 

■ J? , , we obtain 

3{uv) ' 



<8 



y- + l 



(8) 



Case A 
Case A1 
Case A2 
Case A3 
Case B 
Case C 




FIG. 6: (Color online) Support for the generalised Rice 
theorem as a meaningful approximation in turbulent 
channel flows with various Reynolds numbers and wall 
actuations (see Table I), (a) B\ as a function of y+, (b) 
Bi as a function of y/8 with no wall- forcings and (c) B\ 
as a function of y/8 with wall- forcings. 
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where C s is given by 



. x 10 



c, = 



kB 2 C 



(9) 



The classical claims [4[ are that k ~ 0.4, C ~ 2 and 
£?2 — 1 in the intermediate range 1 <C y+ <C 7?.e r as 
lZe T — > oo. These claims therefore imply that C s should 
also be a constant in that same range and limit provided 
B\ is. Whilst, as we have seen, B\ is not too far from 
being constant in the range 8 U <C y < S, k and C are 
significantly far from constant in this range (see Fig. [3] 
and Fig. [7]). Even so, our DNS evidence (see Fig. [5]) 
suggests that C s tends to a well-defined constant in the 
range 8 V y < 6 as 7?.e r increases. Remarkably, this 
condition on 7?.e r for the constancies of C s and Bi seems 
to require as little as lZe T exceeding a few hundred. It 
is equally remarkable that our calculation of N s , which 
underpins B\ and C s , has involved an average over a 
number of time-samples that is two orders of magnitude 
smaller than for the time average required to statistically 
converge (u), (uv), E and e. 




400 



FIG. 7: (Color online) C as a function of y+ for 
different Reynolds numbers and different wall-actuation 
cases (see Table I). 

The constancy of C s in the range b v <C y < 5 implies 
that, in that range, the number density n s of stagnation 
points decreases with distance from the wall as y+ ■ This 
is in qualitative agreement with Fig. [5] which shows the 
stagnation points to be increasingly denser as the wall is 
aproached. 

The constant C s can be interpreted as representing the 
number of turbulent velocity stagnation points within a 
cube of side-length equal to a few multiples of 5„ (sec 
([8])) placed where y equals a few multiples of 10<5„ as 
seen in Fig. (JSh,). This is the lower end of the range 
where n s ~ y^ 1 and seems to be where the upper edge 
of the buffer layer is usually claimed to lie 

Equations ([S]) and @ have been derived by assuming 
well-defined constant values of k, B 2 , C and B\. How- 
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FIG. 8: (Color online) Normalised number of turbulent 

velocity stagnation points for different Reynolds 
numbers and different wall-actuation cases (see Table 
I), (a) C s versus y + \ (b) C s versus y/S. 



ever, our DNS results show that, at the Reynolds num- 
bers considered, B\ and C s are indeed constant but k, 
i?2 and C are clearly not. Equations ((6]) and ([SJ with 
constant dimensionless values of B\ and C s seem to be 
more broadly valid than the assumptions under which we 
derived equation ((SJ. Therefore, in the next section we 
explore the phenomenology behind the new equations © 
and © and the constant values of B\ and C s in the range 
&v *C y ^ S and in section IVTl we go one step further and 
derive the consequences of the constancies of B\ and C s 
on the mean flow profile without assuming well-defined 
constant values of K, B 2 and C. 



V. PHENOMENOLOGY 

One interpretation of the constancy of B\ can be ob- 
tained by considering the eddy turnover time r which 
is defined by e = E/t. Combined with the equation 
e = f which defines A, one obtains 3A 2 = 2vt. Using 
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([6]). B\ = Const is then equivalent to 



-S 2 -£' + — and account should be taken of the fact that 



(10) 



which indicates that in the equilibrium layer, the time 
it takes for viscous diffusion to spread over neighbour- 
ing stagnation points is the same proportion of the eddy 
turnover time at all locations and all Reynolds numbers. 
In high Reynolds number turbulence, the turnover time 
is also the time it takes for the energy to cascade to the 
smallest scales. 

For an interpretation of the constancy of C s note 

first that © and l s = W^vF = K^) -172 imply 

i\ = C~ x 8 v y. From e = § ^ and B 1 = X/i s it then 
follows that 



with 



_ 2 Eu T 

3 K s y 



Bl 
C s ' 



(11) 



(12) 



The meaning of C s and B\ constant is therefore, using 
pip , that the eddy turnover time r = E/eis proportional 
to y/u T throughout the range where they are constant. 
The constant of proportionality is 3k s /2 where n s is de- 
termined by the stagnation point coefficients B\ and C s 
and is constant if they are constant. We refer to t\> s as 
the stagnation point von Karman coefficient. 

Note that, in the present context, equation (fTTj) re- 
places the usual e = u^/ny [i[ , and that these two equa- 
tions reduce to the same one only if and where E ~ u 2 
independently of y+ and lZe r . 



VI. THE MEAN FLOW PROFILE IN THE 
EQUILIBRIUM LAYER 

In this section we spell out the consequences of the 
constancies of B\ and C s on the mean flow profile. 
In the equilibrium layer the expectation is that B 2 — > 1 



in the limit 7Ze T — > oo. This means that 



dy 



B 2 e may be replaced by —{uv)-^(u) = e in the equilib- 
rium layer. The constancy of B\ and C s in this same limit 

3 

implies a constant n s = B\jC s in e = ^E + -^j where 
E + = E/u 2 T . It then follows that -(to)^(u) = e = 

3 

In turbulent channel/pipe flows where we can 
have some mathematical confidence that, as lZe T — > oo, 
— (uv) — > u 2 . in an intermediate layer S v <C y <C S, it 
finally follows that 



dy 



3 K s y 



(13) 



in that same layer and limit. At finite Reynolds numbers 
this new equation (f]~3|) should be replaced by — 



2 . 

3 r K s y 

B 2 , -B3 and k s all have their own, potentially different, 
rates of convergence towards their high Reynolds number 
asymptotic constant values. 

An important step taken in deriving both (d} and (TT3")) 
has been the local high Reynolds number balance V ~ e 
in the equilibrium layer. In terms of the classical assump- 
tion e ~ vi^/ny, V = B 2 e implies that Vy/u^. should 
equal B 2 /k which should be constant in the equilibrium 
layer as a result of the balance between V and e. In 
terms of the new formula (fTTj) . V = B 2 e implies that 
^Py/(E + u^.) should equal B 2 /k s and the balance be- 
tween V and e means that it should be B 2 /k s rather 
than B 2 /k which is constant in the equilibrium layer. 
The main difference is the presence of E in (TTT]) . In Fig. 
[9]we plot our DNS results for ^Py/(E + u^.) versus y+ and 
in Fig. rjD] our DNS results for Vyjv? T . It is clear that 
the collapse between the different Reynolds number and 
wall-actuation data is far worse and the y-dependence 
in the equilibrium layer far stronger for B 2 /n than for 
B 2 /k s . 

These DNS results arc for Reynolds numbers which 
arc not very large; yet the high-Reynolds number con- 
stancy of B 2 /k s in the equilibrium layer seems already 
not exceedingly far from being reached whereas no such 
indication is shown in the plot of B 2 /k versus y+. Fig. 
[TTlis a linear-linear replot of Fig. [9] for easier comparison 
with Fig. H 



-^*— Case A 
-^•^Case A1 
— e— Case A2 
— s— Case A3 
-^Case B 
Case C 




FIG. 9: (Color online) Linear- log plots of versus 
for different Reynolds numbers and different 
wall-actuation cases (see Table I). This is the same as 
B 2 /k s versus y+ because of (fTTj) and V = B 2 e. 

From (|T3"|) . direct plots of f e+u ~d% ( u ) sn0lu d give 1/k s 
in the equilibrium layer when lZe T — > 00 and B 2 /(B^k s ) 
in that layer at finite Reynolds numbers. These plots for 
each of the Table I cases are given in Fig. [T^J They do 
not compare favourably with the plots of ^ {u} , effec- 
tively plots of B 2 /(B^k), in Fig. [3J However, this does 
not mean that in the limit lZe T — > 00, Q is better than 
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Case B 
Case C 



Q. 




FIG. 10: (Color online) Linear-log plots of Vy/u^. 
versus y + for different Reynolds numbers and different 
wall-actuation cases (see Table I). This is the same as 
B 2 /k versus y + because of © and V = B 2 e. 
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FIG. 11: (Color online) Linear-linear plots of B 2 j ' k s 
versus y + for different Reynolds numbers and different 
wall-actuation cases (see Table I). 



(fT3"| in the equilibrium layer. The facts that C s and B\ 
are approximately constant in the range 8 U <C y < 5 and 
that B 2 /n is much less collapsed and less constant along 
y than B 2 /k s at Table I's Reynolds numbers (see Fig. [13] 
and compare it with Fig. [TTj) suggest that the strong y 
and lZe T dependencies of B3 partly cancel those of B 2 j ' k 
at those Reynolds numbers. As the Reynolds number 
is increased to the point where B3 reaches its asymp- 
totic value 1 then this cancellation will either disappear 
if B 2 /k does not tend to a constant or will remain if it 
does. In the specific context of the present stagnation 
point approach, the choice between these two scenarios 
will depend on the high-Reynolds number scalings of the 
kinetic energy E. 

According to classical similarity scalings, as TZe T — > 00, 



- Case A 
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Case A2 
Case A3 
Case B 
Case C 




400 



FIG. 12: (Color online) Plots of versus y+ 

for different Reynolds numbers and different wall 
actuations (see Table I). 
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400 



FIG. 13: (Color online) Linear-linear plots of 
B^^--^j{u) versus y + for different Reynolds numbers 
and different wall actuations (see Table I). These are 
effectively plots of B 2 /k to be compared with the 
similarly plotted B 2 /n s in Fig. [TTJ 



E ~ u^. independently of y, S and v in the equilibrium 
range <5„ <C y <C S. If this is true, then ((4]) and the log- 
law are recovered from (|11[) and (|13[) with a von Karman 
coefficient k proportional to k s . However, Townscnd's 
[23( | idea of inactive motions would suggest that E docs 
not scale as u^. in the equilibrium layer when lZe T — > 00. 
If this is the case, then (fT3"]) does not yield (01 and B 2 /k 
does not tend to a constant in the high Reynolds number 
limit. 

This discussion naturally brings us to the non- 
universality of measured von Karman coefficients Q 
which we now comment on before moving to the anal- 
ysis of some of the highest Reynolds number DNS data 
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currently available. If E ~ u 2 . at high Reynolds num- 
bers and the log-law therefore holds as a consequence of 
(fl~3"|l . then, because of (|12|) . the von Karman coefficient 
will have to be proportional to B\ and inversely propor- 
tional to C s , the number of stagnation points within a 
volume 5f, at the upper edge of the buffer layer. There 
is no a priori reason to expect B\ and C s to be the same 
in turbulent channel and pipe flows, for example. There 
is therefore no a priori reason for the von Karman coef- 
ficient to be the same in different such flows either. 

In the case where the log-law does not hold because 
of the effect that inactive motions have on E + in (fT3"l) . 
data fitted by a log-law may yield different von Karman 
coefficients both as a result of k s = B\ jC s but also as a 
result of fitting mismatches. 

In conclusion, whatever the scalings of E + , one can 
expect measured values of the von Karman coefficient to 
be non- universal as has indeed been recently reported Q . 



VII. HIGH REYNOLDS NUMBER DNS DATA 

We now test some of our results and conclusions on a 
set of data which includes the highest Reynolds number 
channel flow computations currently available [loj . i.e. 
lZe T = 2000. This set also includes data for lZe T = 950 
[Tol | and our own highest Reynolds number DNS data 



{TZe T = 395). We plot 
Fig. HI and JL£{u) 

2^ 3 E + u T 



Ie^tM^ = B ^/(B 3 n s ) (see 
B2HB3K) (see Fig. [T5|) as well as 



B 2 /k s (sec Fig. USD and B 3 ^-±(u) 



B 2 /k (see Fig. Q2). 

These high Reynolds number results support and ex- 
tend the claims made in the previous section: B2/ k s ap- 
pears to have the least departures from constancy in the 
intermediate range, better that B2/(B 3 k) which is how- 
ever better than B2/ (B 3 k s ). The variations of B2/ 'k are 
offset by those of B3 (see Fig. \TE\i which explains why 
B2HB3K) looks better than B2/ \B 3 k s ). The situation 
remains therefore identical to the one we encountered 
with the lower Reynolds number simulations in the pre- 
vious section. The Reynolds number needs to be very 
much larger than 2000 for B3 to come close to its asymp- 
totic value 1 in an intermediate layer, as already shown 
by experimental measurements spanning an ever wider 
Reynolds number range in Q. 

Alternative forms for the mean flow profile at high 
Reynolds numbers have been proposed in the literature 
and we test in Fig. [T^the suggestion of a power-law form 
[13, HH . The high Reynolds number data we are using 
here appear to give significant support to such a power- 



law form with power exponent n 



V d j j 

U+ dy U + 



2/15, 



y + ' 1+ "^ 15 ^ in the intermediate layer. On 



i.e. -F-U+ 

dy + 

the basis of (fT3"|) . this result suggests that E + has a 
power-law dependence on y + in that same layer. In- 
deed, combining (|1 3|) in its finite Reynolds number form, 
§|2£ w ith d U+ 

3 B 3 t Ks y * dy + 




FIG. 14: (Color online) 



» 



as function 



2 E+u T dy \ ' B3K S 

of (a) y+ and (b) y/S. DNS of turbulent channel flows 
without wall actuations. The lZe T = 950 and 2000 data 
are from fioj 



yields 

2/15 



\Ba, i.e. a constant value of 



*<«> 



B 4 ,,-(1+2/15) 



E + y^_ ^ B2 / B% in the equilibrium layer if S4 is constant 
in that layer. Figure [20] supports this conclusion albeit 
with a constant value of E+y 2 / 15 B 2 / B3 which appears 
to increase slowly with Reynolds number. This Reynolds 
number dependence may be intrinsic to E + resulting, 
perhaps, from Townsend's inactive eddy hypothesis. 



VIII. CONCLUSION: A NEW ASYMPTOTIC 
APPROACH 

Our DNS suggest that B\ = X/£ s and C s = n s 8^y+ are 
approximately constants in the region 5 V y < 5. Our 
DNS demonstrate these well-defined approximate con- 
stancies for lZe T as low as a few hundred. These con- 
stancies imply that, in the region 8 V <C y ^ 8, the eddy 
turnover time r equals ^K s y/u T with k s = B\jC s . As- 
suming the constancies of B\ and C s to be early man- 
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FIG. 15: (Color online) ^rg^(w) = ^ as function of FIG. 16: (Color online) lB 3 ^\-^(u) = as 

(a) y + and (b) y/S. DNS of turbulent channel flows function of (a) y + and (b) y/S. DNS of turbulent 

without wall actuations. The lZe T — 950 and 2000 data channel flows without wall actuations. The lZe T = 950 

are from [lj|; these plots have have already been and 2000 data are from [l(| 
presented in [l(| for this data 



ifestations of a high Reynolds number behaviour, i.e. 
that B\ and C s and therefore k s are independent of 
y in 6 U <C y < d as TZe T — > 00, it then follows that 
= ^E+jr^ in the equilibrium region 8 U <C y <C 5 
where production may be expected to balance dissipation 
and —{uv) ~ u\. The asymptotic equality — (uv) ~ v% 
is mathematically supported only for turbulent chan- 
nel/pipe flows. 

The classical intermediate asymptotics approach, 
which assumes no dependence of the mean velocity gra- 
dient on v and 6 where b v <C y <C 5, does not consider 
the effect that Townsend's inactive motions [23| may or 
may not have on -^(u). However, if a new intermediate 
asymptotic approach is taken where the assumption of 
no dependence on v and S is applied to r = E/e instead 
of ^(w), we are then led to r ~ y/u T: and the small ef- 
fects of inactive motions on ^- (u) may not be neglected 



because r ~ ^n s y/u T . —(uv) ~ v? T and -(jto)^(ii) ~ e 
yield (fTU)) which explicitly contains E + . If as a result 
of inactive motions, E does not scale as u\ then this re- 
vised intermediate asymptotic approach will not predict 
a log law for the mean profile even though a stagnation 
point von Karman coefficient k s exists and is well de- 
fined within the approach. The mean flow prediction of 
this approach is instead controlled by the intermediate 
asymptotic dependence of E + on y + and lZe T . If this 
dependence on y + is a power law E + ~ y7 n in the inter- 
mediate range <5„ <C y <C 5, then the mean flow profile 
will also be a power law, i.e. -^-U^ ~ V+ ln , in that 
intermediate layer. 

DNS of turbulent channel flow with the highest values 
of 1Ze T currently available [To[ suggest the same n = 2/15 
in both E + ~ y7_ n and ^~U + ~ in support of 

our procedure and formula (fT5)) . However, we should 
caution against extrapolating this value of the exponent 
n to higher values of Re T , in particular in the laboratory 
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FIG. 17: (Color online) B 3 ^£^(u) = ^ as function of 

(a) y + and (b) y/8. DNS of turbulent channel flows 
without wall actuations. The lZe T = 950 and 2000 data 
are from [l(J 



A 0.6 




FIG. 18: (Color online) B 3 = —{uv)/u^ as function of 

(a) y + and (b) y/S. DNS of turbulent channel flows 
without wall actuations. The Tie T = 950 and 2000 data 
are from flOl 



where the boundary conditions are in fact different from 
the DNS which is periodic in x and z. We stress that 
the main point of value here for us is the support that 
these exceptional DNS [l(| bring to our approach and in 
particular to our formula (|13p . 

Note also that the stagnation point von Karman coef- 
ficient is defined by e ~ \E+^ J - in the range 8 V <C y < 8 

o k s y 

irrespective of whether the mean flow profile is a log law 
or a power law. Power law profiles of E+ and 
in the intermediate layer 8 V <C y <C 8, simply force 
e ~ %E + ^ J - to imply that the classical relation e ~ — 
does not hold in that layer. 

This new intermediate asymptotic approach is sup- 
ported by our DNS observations that ((6|) and dHJ are 
valid in the region 8 V <C y < 5 because © and © 
imply r ~ ^K s y/u T with k s = Bf/C s . The relation 
K s = Bf/C s offers a link between the underlying flow 
structure, described in terms of stagnation points, and 



the dissipation/cascade statistics of the turbulence. The 
universality in terms of both lZe T and flow-type depen- 
dencies of k s becomes a question concerning the univer- 
sality of the stagnation point structure of the turbulent 
fluctuations. To what extent does it depend on boundary 
and wall forcing conditions? Is it the same in turbulent 
channel and turbulent pipe flows? Is it the same in DNS 
of such flows where periodic boundary conditions are 
used and in laboratory realisations of such flows where 
boundary conditions are clearly not periodic? These are 
questions which must be left for future investigation, but 
our approach makes them fully legitimate as there is no 
reason to expect the stagnation point structure of tur- 
bulent fluctuations to be exactly the same in all these 
cases. 

The implications of this new approch for the mean flow 
profile in turbulent channel/pipe flows come by invoking 
a local balance between production and dissipation as 
well as — (uii) ~ u\ in the intermediate range 8 V <C y <C 8 
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FIG. 19: (Color online) Power law mean velocity profile: 

n = ij^-^U+ plotted against y + . DNS of turbulent 
channel flows without wall actuations. The lZe T = 950 
and 2000 data are from [lj| 



as lZe 7 
2 e 



oo. A direct test against data of ^j— U+ 



+ in that same range and limit cannot be expected 
to be successful if lZe T is not large enough for —(uv) to 
equal u 2 T over the range 5 V <C y <C S. As clearly shown 
by various experimental and numerical data, this equal- 
ity is well beyond the highest Reynolds numbers cur- 
rently available both numerically and in the laboratory. 
The siginificant finite Reynolds number deviations from 
— (uv) ~ u\ appears to compensate the deviations from a 
log-law and from the local production-dissipation balance 
with the result that plots of y-^U+ have a less varying 



appearance than plots of 



At face value this 



could be misinterpreted as better support for the log-law 



_d 



— than for 

dy + 



2 



3 K s y + 



However, plots 



of Bij ' k s look significantly better than plots of B 2 /n thus 
demonstrating that f^g^C^f = B 2 /(B 3 n s ) looks worse 

than y-^jU + = B 2 /{B 3 k) only because B 3 = —(uv)/u* 
is so significantly non-constant (see Fig. fT8)) . 
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FIG. 20: (Color online) (a) E+ and (b) .E+y+ff with 
n = versus y+ for DNS of turbulent channel flows 
without wall actuations. The lZe T = 950 and 2000 are 
from [ill 
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As a final note, future investigations should attempt to 
uncover the small-scale intermittency corrections to our 
scalings by determining the weak dependencies that B\ 
may have on lZe T and y + as a result of small-scale inter- 
mittency. These dependencies will cause dependencies of 
k s on TZe T and y+. It is remarkable that small-scale in- 
termittency may have an impact, even if small, on the 
scalings of mean flow profiles. These Reynolds number 
parts of their scalings should be distinguished from those 
that Townsend's attached eddies may be contributing via 
the scaling of E + on lZe T . 



Appendix: Numerical method for the computation 
of stagnation points 

In this study of turbulent channel flow, we focus on 
stagnation points of the fluctuating velocity field, i.e. 

u'(x, t) = u(x, t) - (u) = (A.l) 

where (.) here, denotes an average in space over the 
homogeneous directions x and z at a particular time. 
These zero- velocity points are Galilean invariant and re- 
sult from the intersections of the three random surfaces 
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u'(x,t) = 0, v'(x,t) — 0, w'(x,t) — 0. Intersections of 
two random surfaces gives lines and the intersections of 
these lines with a third random surface gives points. 

A root finding method is required to obtain where in 
space our random function u'(x,t) is locally zero. Here 
we use the iterative Newton-Raphson method, 

Xnew = x Q id + Sx with [Vtt'] £ Sx = —u' c (A. 2) 

where Vti'(x, t) = Vu(x, t) — (Vu) and the subscript £ 



stands for an interpolated quantity. The interpolation of 
the velocity and its gradient was done using fourth-order 
Lagrangian interpolation [26j]. The particular choice of 
the interpolation was chosen based on robustness, accu- 
racy and computational efficiency. The 3x3 linear system 
of equations was solved for Sx simply using the Cramer's 
rule 



Sx 



det(«£, d Q y \ c , S q z | £ ),det( \ c ,u' c , 9 g z | £ ),det( g" \ c , °™ y \ £ ,u' c ) 



det(^| £ , 



du' 

dy 



du' 



(A.3) 



assuming that det(^| £ , %\ c , ^ 0. 

The Newton-Raphson method gives a very rapid local 
convergence to a root if the initial guess is sufficiently 
good. So to have a chance at good starting points we 
took them midway between two consecutive grid points 
throughout our computational domain. It is well known 
that different initial guesses can converge to the same so- 
lution, due to the unpredictable global convergence prop- 
erties of this iterative method. To avoid this issue we try 
to bound the root finding no further than the neighbour- 
ing computational cells and we also ensure that no more 
than one stagnation point exists in a cell, which is what 
should be expected from a smooth velocity field of a well 
resolved DNS. 

Our method was also compared to an algorithm pro- 



posed in [27[ for the detection of unstable periodic or- 
bits in chaotic dynamical systems, which has good global 
convergence due to its attracting nature. This method is 
based on a universal set of linear transformations, which 
transform unstable fixed points to stable ones whilst 
mainting their positions. However, this method can be 
expensive in more than two dimensions and this is the 
reason why we choose here Newton-Raphson, which is as 
accurate. 

The number of zero-velocity points N s are computed 
within thin slabs of dimensions L x x S y x L z , with 5 y ~ S v , 
parallel to the channel's wall. Time averages of N s were 
taken by repeating the same procedure for several time 
instances. 
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